Pre-melting of crossing vortex lattices 
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The pre-melting of high vortex density planes observed recently in layered superconductors in 
tilted magnetic field l] is explained theoretically. Based on the structural information of the crossing 
lattices of pancake and Josephson vortices the effective vortex cage potential at different lattice 
sites is determined numerically. Melting takes place when the thermal energy allows proliferation 
of vacancy-interstitial pairs. It is found that the increased density of pancake vortex stacks in the 
planes containing Josephson vortices, rather than their incommensurate structure, is the main cause 
for pre-melting. 
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Melting of solids, an everyday phenomenon from the 
shrinking ice cube in a drink to the processes in the earth 
mantle, is far from being well understood. At the melt- 
ing transition the shear modulus of the material vanishes 
[2]. There are several phenomenological criteria predict- 
ing when this will happen (see [2H4] for a discussion). 
Most prominent is the Lindemann criterion [5], stating 
that melting occurs when the mean square fluctuation 
of the atoms in a solid is a fraction cl (the Lindemann 
number) of the interatomic distance. On a more concep- 
tual level it is expected that the proliferation of defects 
or grain boundaries plays an important role in melting. 
In two dimensions melting was indeed shown to be a two 
stage process, driven by unbinding of dislocations and 
disclinations, respectively [6], while in three dimensions 
a sudden increase of the vacancy concentration was ob- 
served at the transition [3j [?H9] . 

In this Letter we address the question of the melting pro- 
cess in a heterogeneous system. Since inhomogeneities 
and extended defects occur naturally in any practical sys- 
tem, comprehension of the melting transition in such sys- 
tems is of broad importance. Vortex lattice in supercon- 
ductors has been extensively utilized as a model system 
for theoretical and experimental studies of the melting 
transition [TO] [TTj . Melting in flux line lattices of type- 
II superconductors is characterized by the dimensionless 
parameter ex = A/A T . A denotes the London penetration 
depth and A T is a thermal length scale, only related to 
the temperature T and the flux quantum (/>o = hc/(2e), 
A T = (j)l/{1^7i 2 k B T). In high-T c materials with their 
elevated transition temperatures and large values of A, 
ex can become of order one replacing a large part of the 
mean field phase diagram by the vortex liquid phase [12] . 

It was recently found pQ that by tilting the magnetic 
field away from the c-axis the melting transition in highly 
anisotropic high-T c superconductor Bi 2 Sr2CaCu208+<5 
(BSCCO) changes from a homogeneous first-order pro- 
cess into a two-step transition with an intermediate 
lamellar solid-liquid phase that is apparently driven by 



an intrinsic heterogeneity of the vortex lattice structure. 
In layered materials like BSCCO the vortex lines along 
the c-axis consist of stacks of pancake vortices (PVs) re- 
siding in the superconducting layers, whereas the vor- 
tex lines parallel to the ab plane take advantage of the 
normal layers forming Josephson vortices ( JVs) [13] [14] . 
When the magnetic field is tilted away from the c-axis, 
chains of higher vortex density were found in Bitter deco- 
rations [T5HT7] , scanning Hall probe microscopy p~8j [19] , 
Lorentz microscopy [20], and magneto-optical imaging 
[2TH23] , and described as crossing lattices of pancake and 
Josephson vortices (see Fig. la) [24 - 30 . The resulting 
vortex lattice structure displays a very rich phase dia- 
gram that depends critically on the strength and direc- 
tion of the magnetic field H and on the ratio of e s /7, 
where e s = X/s (e s ~ 130 in BSCCO), 7 = A c /A is the 
anisotropy parameter (7 ~ 500), A = A a 5 is the in-plane 
penetration depth, and s is the interlayer spacing [29J. 
The essential element that characterizes the crossing lat- 
tices state is that the PV stacks that reside along the JV 
planes (see Fig. la) experience local environment that is 
substantially different from the bulk PVs thus leading to 
formation of a heterogeneous system which is expected 
to undergo a distinctive melting process. Since the de- 
gree of heterogeneity can be readily controlled by the 
direction and magnitude of the magnetic field, the study 
of this system could provide an important insight into 
the microscopic mechanisms of melting in heterogeneous 
systems. 

The previous theoretical studies of the melting of the 
crossing lattices [24-29 and most of the experimental 
studies [3TU33] have investigated the global behavior as- 
suming a single melting transition. The recent magneto- 
optical imaging of the melting process has shown, how- 
ever, that the local process in such a heterogeneous sys- 
tem cannot be described by a single transition and that 
the planes containing the Josephson vortices melt first 
creating an intermediate periodic solid-liquid state pQ. 
The goal of the present work is to provide a theoreti- 
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cal framework describing the melting of a heterogeneous 
lattice and of the pre-melting of the Josephson planes. 
Keeping in mind the complicated lattice structure as well 
as the absence of a comprehensive theory of 3D melting 
even for simple systems, we will resort to a mean- field-like 
vacancy theory of melting. In this description the effec- 
tive cage potential, acting on a given PV due to its inter- 
action with all other vortices, is calculated by replacing 
the instantaneous positions of the latter by their equilib- 
rium ones. Melting is concluded to occur at temperatures 
where the PV under consideration can escape from the 
cage potential via its lowest saddle point, triggering pro- 
liferation of vacancy-interstitial pairs. This happens first 
for PVs in the Josephson planes. The melting temper- 
ature determined in this way is in good agreement with 
the experimental findings. 

Model. — We follow here the conclusions of [29 that the 
discrete layer structure has strongest influence on the 
cores of tilted and Josephson vortices, but that interac- 
tion contributions to the total energy usually can be com- 
puted within continuous approximation. Pancake and 
JVs in their equilibrium position can be considered to be 
threaded on a contour line r^(t). The vorticity uj of the 
lattice can be written as 

w(r) = j dt{dvi{t)/dt)5 (r - n(t)) , (1) 

i 

which is related to the magnetic induction B by 

B + Vx^ A^(V x B) a f a = -0o"(r). (2) 

a 

Here X z = A c and X x = X y = A. The free energy T of 
vortex systems is a functional of the vorticity field which 
includes contributions both from the magnetic and the 
Josephson interaction [29] 

4>l f d 3 k K| 2 + fc 2 (AgK k | 2 + A 2 |a, ± k| 2 ) 
M 8irJ (2tt) 3 (l + A 2 fc 2 )(l + A 2 fc 2 +A2k 2 L ) ' 

(3) 

u?k is the Fourier transform of u;(r). The subscript _L 
denotes the projection of a vector onto the x?/-plane. The 
actual equilibrium flux lattice structure follows from the 
minimization of 

- J d 3 rB • H/4tt (4) 

for a given magnetic field H. Its determination is in gen- 
eral a complicated task because of the potentially large 
number of competing ground states. In this paper we 
will not address this question but instead assume that 
the structures found experimentally are indeed the free 
energy minima for a given applied magnetic field H. Fur- 
ther below we will consider the effects of a possible incom- 
mensurate superstructure and find that these are weak. 




FIG. 1. (a) Schematic 3D plot of the crossing PV-JV lat- 
tice. The JVs are aligned along the x axis and are stacked 
along z forming "Josephson planes" along xz-planes. (b) Top 
view (x2/-plane) of the PVs (red) and JVs (blue). The vortex 
numbers are used in Figs. 2 and 3. (c) Front view of the 
JV lattice (^/z-plane). (d) Side view of PVs and JVs in the 
Josephson plane (xz-plane). 



Once the equilibrium structure is known, the effective 
potential for a pancake vortex displacement u^ n in stack 
i and layer n can be found from the change of the vor- 
ticity field Au?(r) in (J3|. Its Fourier transform is Au?k = 
zse lV - Vi [e zk ^ Ui '™ — l] . This allows to calculate the free 
energy change AJ r (n^ n ) which depends in general on 
i and n. The actual vortex structure in the Joseph- 
son planes consists of continuous vortex lines formed of 
stacks of PVs, connected by JVs, as indicated in Fig. [T]i 
[29], 134] . The resulting discontinuity of the JVs in a sin- 
gle plane (and hence also the displacements of the PVs 
from a straight stack) is of the order A 2 /[(2n — 1)75] [34], 
which is for our parameter values < 0.06 a j/n. n > 1 
denotes the number of layers between the JVs [34 i. We 
will therefore ignore in the following this small effect and 
assume a simple crossing lattice of straight stacks of PVs 
and straight JVs [35] . 

Displacing a single PV from a straight stack gener- 
ates an antiparallel pair of JVs in the adjacent non- 
superconducting layers. Because of their large mutual 
annihilation, the interaction of this excited JV pair with 
the static JVs of the crossing lattice will be very small, 
provided the displaced PV is not too close to the JV. We 
will ignore this weak interaction in the following as well. 
With these approximations Co of PV and JVs are now 
perpendicular to each other, and hence Josephson vor- 
tices do not contribute to the effective potential of the 
displaced pancake vortices once the structure is fixed. 
Moreover, the effective potential for a single displaced 
PV does not depend on its position in the stack. This 
approximation is justified the more the larger the dis- 
tance between the displaced PV and the JV. 
After summation over the PVs in the stacks and sup- 
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FIG. 2. Set of contour plots of the cage potential AT(ui) for 
a single displaced PV in a stack of PVs in the Josephson plane 
for the different stacks numbered as in Fig. lb. Only the PV 
with rest position in the center is assumed to be displaced, 
keeping all other PVs at their equilibrium positions. With 
the choice of aj/a = 5/6 this sequence is repeated every 6 
stacks. 

pressing the layer index n, the effective potential AJ r [n i ] 
for a single PV, displaced from a straight stack by i^, 
can be written in the form 



e T e t 



[K (\K ij -u i \/X)-K (R ij /X)}. (5) 



Here Kq(x) is the modified Bessel function, R^- the po- 



sition vector connecting stacks i and j, and Ri 
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The remaining sum has to be performed numerically, 
which requires the knowledge of the flux lattice struc- 
ture. We will make some assumptions for the ground 
state, concluded both from analytical theory and exper- 
imental findings [15J [29] : 

(i) All vortices are located in planes parallel to the xz- 
plane, assuming H y = (Fig. 1). Planes including cross- 
ing PVs and JVs are called Josephson planes. 

(ii) PV stacks outside the Josephson planes are assumed 
to form a equilateral triangular lattice of spacing a p~6j 

Eg. 

(iii) The PV distance in the Josephson planes is denoted 
by aj < a and the distance between Josephson and the 
adjacent planes is denoted by b > \^3a/2 (Fig. lb), aj < 
a is the result of the attractive interaction between PVs 
and JVs. 

(iv) The projection of the JVs onto the ?/z-plane forms 
squeezed Abrikosov lattice with the lattice parameter ra- 
tio bj/c = V3j/2 (Fig. lc). 

The above conditions can be somewhat relaxed as dis- 
cussed below. 

To be specific we used the lattice structure found in the 
Bitter decoration experiment [15 performed in BSSCO 
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FIG. 3. Contour plot of A J 7 for PVs in the plane adjacent to 
the Josephson plane (numbered 7 to 11 in Fig. lb) and in the 
bulk (number 12). Melting occurs at about 10 K higher than 
for PVs in the Josephson plane in Fig. 2. 



at H = (30,0, 24) Oe where a = 1, aj = 0.833 , b = 
0.91, c = 0.04, bj = 17.32, A(T) = A(0) y/\ - T/T c , A = 
0.2, all lengths in fim. The value of H x is slightly larger 
than that used in [1 j. Figure [2] shows the results of the 
calculation of the free energy change AJ 7 due to the dis- 
placement of the individual PVs in the Josephson plane 
labeled 1 to 6 in Fig. lb. Since in the present case 
a j j a = 5/6, the potential shapes recur every six vor- 
tices. The contour maps of lines of constant energy are 
depicted in units of ksT. There are two important ob- 
servations here. First, although the precise shape of the 
potential is different for the individual vortices, saddle 
points with very similar value T « 75 K appear for 
all the locations. At this temperature PVs can escape 
in direction transverse to the Josephson plane, vacancy- 
interstitial pairs proliferate, and the lattice starts to melt 
along the Josephson planes. This transverse melting pro- 
cess appears to be uniform for the various PVs along the 
Josephson plane. The second observation is that in con- 
trast to the transverse direction, the shape of the poten- 
tial along the Josephson plane is rather position depen- 
dent. Vortices that are located at symmetric points that 
are in registry with the adjacent planes, like vortex 3 and 
6, show potential that is quite localized in the x direc- 
tion. On the other hand, the potentials at asymmetric 
locations, like vortex 1 and 5, show extended protrusions 
that indicate enhanced vortex mobility at the less-stable 
vortex positions along the Josephson planes. This finding 
can explain the apparent observation of site dependent 
enhancement of vortex fluctuations along the Josephson 
planes in Lorentz microscopy studies [20 that was dis- 
cussed in terms of the incommensurability of the vortex 
chains structure. 
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The corresponding pictures for PVs in the plane adja- 
cent to the Josephson plane are shown in Fig. 3 along 
with PV number 12 in the bulk. Despite of the fact that 
each of the PVs has a different local energy landscape, it 
is clearly seen that all the PVs can leave their positions 
only at temperatures of T m c± 85 K, 10 K higher than 
the PVs in the Josephson planes. The experimental bulk 
melting temperature is T m « 90 K. A remark is in order: 
since the vortex cage potential is calculated with all sur- 
rounding vortices in their equilibrium positions, prolifer- 
ation is predominant in the direction perpendicular to the 
Josephson plane. We expect that in a treatment which 
allows simultaneous motion of all vortices this prevalence 
is reduced. 

Figure [4] shows a detailed calculation of the melting tem- 
perature T m of the PVs in the Josephson planes as a 
function of their intervortex distance aj and the separa- 
tion b between the Josephson plane and the adjacent PV 
plane (see Fig. lb). T m is found to drop rapidly with 
increasing a/aj and with decreasing b/a giving rise to 
pre-melting of the Josephson planes that could in princi- 
ple be as large as 20 K. In reality, however, aj and b are 
not independent parameters and are both determined by 
H x and H z . In particular, increasing H x decreases aj 
but increases 6, so that the two dependences moderate 
each other to a large extent. We therefore expect a much 
smaller pre-melting as observed experimentally pQ. 
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FIG. 4. Melting temperature of PVs along the Josephson 
planes T m in crossing lattices state with regular lattice pa- 
rameter a = lfim as a function of b/a (left panel) and of a/ a j 
(right panel). 



So far we assumed equidistant vortex positions in the 
Josephson plane. This assumption can be somewhat re- 
laxed by allowing a modulation of the PV spacing, i.e. 
a j — » aj(n) where n numbers the intervals. The aj(n) 
follow from the competition between the PV-JV interac- 
tion in the Josephson plane on one side, which favors a 
lattice spacing smaller than a, and their interaction with 
the PVs outside the Josephson plane on the other side, 
which favor registry with the lattice constant a, which 
we assume to be fixed. Such a system can be described 
by the Frenkel-Kontorowa model [36 , 



H 



= £[(« 



„+l -0 n -6) -2( COS r , 



(6) 



The n describe the modulation of the PV positions, 



x n = an+a6 n /(27r), £ = 2tt 2 &2 / (q 2 ki) , and q = a/A. The 
misfit parameter 5 < 0, depending on H x , favors a higher 
PV density in the Josephson planes. It will be determined 
below from the average vortex spacing (aj(n)} = aj. 
The coefficient k\ follows from the vortex interaction as 
ki = Ko (q) + q~ 1 Ki (q) . k<i is determined numerically by 
fitting the actual vortex interaction of a Josephson plane 
stack with its neighboring stacks in the satellite plane to 
a cosine model, taking up to the fifth next neighbour into 
account. The ground state configuration follows from the 
saddle point condition 9 n +i + — 2# n = (sm9 n . In 
the continuum limit, n — » 0(n), one arrives at the rigid 
pendulum equation 0"{n) = (sin 6. Its solution can be 
expressed by elliptic functions which depend on the con- 
stant of integration rj. r] follows from the minimization 
of K(r] 2 - 2) + 4E - TRfy/VC Here K(rj) and E(rj) are 
the complete elliptic integrals of first and second kind, 
respectively. For small misfit, \S\ < S c = A^/^/tt, 0(n) 
locks-in at a multiple of 27r, corresponding to aj = a. 
For \S\ > \S C \ the solution for 6{n) is staircase-like with 
horizontal terraces at 0(n) ~ 2np (p G N), which are 
connected by steps of width ~ C -1 ^ 2 - The average PV 
spacing in the Josephson plane is in general incommen- 
surate with a. For a given aj, 5 is determined from the 
relations aj = a [l - y/C/(2r)K)] and rj = 4^E/(ttS). 
With the parameter values given above we get £ = 0.409, 
\5 C \ = 0.814, r] = -0.876 and S w -1.115. Thus the sys- 
tem is indeed in the incommensurate phase . The result- 
ing modulation of the vortex distance aj(n), however, 
is less than 6% of the mean vortex distance and hence 
can be safely ignored for the investigation of the melting 
transition. 



In summary, using a cage potential model we have cal- 
culated the melting temperature of the PVs at differ- 
ent locations in the crossing lattices state parametrized 
by three different lattice constants a,aj,6. Detailed re- 
sults were obtained for a given PV lattice structure rel- 
evant to highly anisotropic superconductor BSCCO. We 
showed that melting of crossing flux line lattices sets in 
along the planes containing JVs where the local melting 
temperature can be substantially lower than in the bulk 
resulting in a periodic solid-liquid structure as observed 
experimentally [1 . We have analyzed the effect of the 
incommensurate modulation of the PV lattice constant 
aj along the Josephson plane and found it to be too 
weak to have a significant effect on the melting transi- 
tion. The primary cause of the pre-melting transition is 
proliferation of vacancy-interstitial PV pairs at reduced 
temperatures due to the increased PV density along the 
Josephson planes. 
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